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1.  Introduction 


Much  of  the  literature  that  presents  analytical  or  computational  solutions  to  “impact-type”  initial 
boundary  value  problems  assumes  an  impact  boundary  condition  in  the  form  of  a  constant  applied 
step  in  stress  or  velocity,  cf. 1  3  This  simplifying  assumption  is  particularly  useful  for 
characterizing  elastic  shocks  in  flyers  and  targets  that  are  homogeneous,  linear  elastic,  and 
semi-infinite.  Here,  elastic  shocks  are  defined  by  stress  wave  amplitudes  that  are  less  than  or 
equal  to  the  Hugoniot  elastic  limit,  cf.,  page  248,  equation  (10.4)  of  Davison.4  However,  for  a  true 
impact  into  a  piezoelectric  medium  that  is  multilayered,  or  inhomogeneous,  or  possesses  other 
material  discontinuities,  the  target  introduces  multiple  reflected  or  scattered  waves  at  the  impact 
face,  so  that  the  history  of  the  stress  a  and  the  particle  velocity  v  at  the  impact  face  is  generally 
time-dependent  and  cannot  be  known  in  advance  of  the  solution  to  the  problem.  The  use  of  a 
constant  stress  or  velocity  boundary  condition  is  an  invalid  assumption  for  modeling  impact  in 
such  targets,  yet  its  use  is  widespread  in  the  mechanics  literature. 

Both  Redwood5  and  Steutzer6  study  the  problem  of  multiply  reflected  transient  waves, 
characterized  by  an  infinite  set  of  jump  discontinuities,  in  linear  piezoelectrics  subjected  to  either 
a  step  in  stress  or  voltage,  but  these  boundary  conditions  are  not  appropriate  for  the  elastodynamic 
impact  problem.  Other  works  address  shock  wave  propagation  in  nonlinear  piezoelectric 
media,7-9  but  the  complexity  of  the  material  response  prohibits  an  analytical  treatment  of  the 
electromechanical  coupling  in  the  unstressed  region  ahead  of  the  shock  wave.  In  this  report,  we 
first  outline  the  linear  equations  of  dynamic  piezoelectricity,  following  closely  the  development 
found  in  Redwood5  and  Steutzer,6  and  show  how  the  d’Alembert  solution  is  also  applicable  to  the 
study  of  multiply  reflected  transient  waves  in  linear  piezoelectrics.  We  then  study  the  “true” 
elastodynamic  impact  problem  using  a  new  impact  boundary  condition  originally  introduced  by 
Scheidler  and  Gazonas;10  the  full  derivation  of  the  impact  boundary  condition  is  based  on  the 
d’Alembert  solution  of  the  wave  equation,  and  appears  in  Gazonas  et  al.11  Our  results  are 
analytical  in  the  Laplace  domain  but  are  numerically  inverted  to  the  time  domain  using  a  modified 
Dubner-Abate-Crump  (DAC)  algorithm12  designed  to  mitigate  the  effects  (using  Lanczos 
a-factors)  of  the  Gibbs  phenomenon  (ringing  seen  in  Fourier  series  representations  of  jump 
discontinuities).  One-dimensional  (1-D)  impact-induced  stress,  velocity,  and  electrical 
displacement  fields  at  the  center  of  a  lead  zirconate  titanate  (PZT-4)  target  are  compared  with 
solutions  obtained  using  a  finite-difference  time-domain  (FDTD)  code  and  the  commercial  finite 
element  (FE)  code,  COMSOL  Multiphysics.13  The  FDTD  code  was  developed  by  Raymond  A 
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Wildman  and  the  COMSOL  simulations  were  performed  by  David  A  Hopkins  to  compare  with 
the  derived  numerical  inverse  Laplace  transform  (ILT)  solutions.  The  final  value  theorem  (FVT) 
can  be  used  to  verify  that  the  asymptotic  (long-time)  stress  and  velocity  in  the  piezoelectric  target 
are  independent  of  the  target’s  elastic  or  piezoelectric  properties,  and  only  the  elastic  properties  of 
the  flyer  and  half-space,  a  result  also  proven  using  the  time-domain  solutions  obtained  for  purely 
elastic  targets.1 1  The  FVT  is  also  used  to  determine  a  closed-form  expression  for  the  asymptotic 
displacement  current. 


2.  Plane  Waves  in  Piezoelectric  Media 


The  1-D  constitutive  equations  for  a  linear  piezoelectric  material  can  be  written  using  the 
notation5, 14 


a(x ,  t)  =  CD 


du{x ,  t ) 
dx 


hD(x,  t )  , 


(1) 


£{x,t)  =  +D(x,t)/e  ,  (2) 

ox 

where  a  is  the  stress  in  N/m2,  CD  is  Young’s  modulus  measured  at  constant  D  in  N/m2,  u  is  the 
particle  displacement  in  m,  D  is  the  electric  flux  density  in  C/m2,  e  =  e0  er  is  the  dielectric 
permittivity  in  F/m,  e0  is  the  permittivity  of  a  vacuum  in  F/rn,  er  is  the  relative  permittivity  of 
the  medium,  h  is  a  piezoelectric  constant  measured  at  constant  D  in  V/rn.  and  £  is  the  electric 
field  in  V/m  (or  N/C ).  The  piezoelectrically  stiffened  elastic  constant  C D ,  Eq.  1,  is  defined  in 
Cady15  (pg.  270  equation  (271  a)),  by  The  Institute  of  Radio  Engineers  (IRE)14  (pg.  1394,  Table 
III)  or  even  Chen  et  al.1  (pg.  4760),  and  is  given  by 

CD  =  CE  +  eh,  (3) 


where  CE  is  Young’s  modulus  measured  at  constant  £,  e  is  a  piezoelectric  coupling  constant  in 
C/m2,  and  h  =  e/e.  In  addition,  from  Newton’s  second  law, 


d2u(x,  t ) 
P  8t2 


d a(x, t ) 
dx 


for  t  >  0,  x  G  (— oo,  oo)  , 


(4) 


where  p  is  the  mass  density  in  kg/m3.  Following  the  development  in  Redwood5  or  Steutzer6  for 
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linear  piezoelectric  media,  substitution  of  Eq.  1  into  Eq.  4  gives 


d2u(x,t)  CD  d2u(x,t)  h  dD(x,t) 

d t2  p  dx2  p  dx 

If  there  is  no  free  charge,  pfree  =  0,  inside  the  piezoelectric  medium,  then  div  D  =  pjree  =  0  (see 
e.g.,  Lorrain  and  Corson,16  pg.  109),  and  the  differentials  with  respect  to  y  and  z  are  assumed  to 
be  zero  (plane  wave  propagation  in  Cartesian  coordinates),  dDy/dx  =  dDz/dx  =  0,  hence 
d,Dx/dx  =  0,  and  Eq.  5  becomes 

d2u(x,t )  2  d2u(x,t ) 

dt 2  dx2  ’ 

where  c  =  \JCD / p  is  the  piezo-modified  longitudinal  wave  speed  with  CD  given  by  Eq.  3. 


3.  Impact-Induced  Plane  Waves  in  a  Piezoelectric  Target 


We  now  consider  the  problem  of  a  semi-infinite  non-piezoelectric  elastic  flyer  traveling  at  initial 
velocity  V0  that  collides  with  (and  welds  to)  a  finite,  undeformed,  stationary  piezoelectric  target  of 
finite  length  /;  the  target  is  backed  by  a  semi-infinite  non-piezoelectric  elastic  half-space  (Fig.  1). 
Because  our  sign  convention  for  stress  is  taken  positive  in  compression,  the  stress  in  the  linear 
elastic  piezoelectric  target  is  given  by  Eq.  1  with  the  minus  sign: 

<?i (x,t)  =  -Ci  dUl +  hDi(x,  t)  .  (7) 

ox 


We  denote  field  quantities  with  a  subscript  of  “1”  as  quantities  in  the  piezoelectric  target,  and 
subscripts  of  “0”  and  “h”  refer  to  field  variables  in  the  flyer  and  half-space,  respectively.  From 
Eq.  6  with  this  notation,  the  transient  target  displacements  Ui(x,  t )  are  governed  by  the  wave 
equation  in  one-spatial  dimension, 


d2ui  (x,  t )  2  d2u\ (x,  t ) 

dt 2  1  dx 2 


for  t  >  0,  x  G  (0,  l )  . 


(8) 


The  general  solution  to  the  wave  Eq.  8  in  the  target  can  be  constructed  using  d’Alembert’s 
equation: 

ui(x,t)  =  fi(t  -  xf ci)  +  gi(t  +  x/ci)  ,  (9) 

where  /i(-)  and  gi(-)  are  arbitrary  functions.  The  displacement  u \  is  assumed  to  be  continuous 
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Fig.  1.  A  semi-infinite  non-piezoelectric  elastic  flyer  (impedance  zq)  traveling  at 
initial  velocity  Vo  impacting  a  piezoelectric  target  (impedance  zi)  of 
length  l:  the  short-circuited  target  V\  (t)  =  0  is  backed  by  a  semi-infinite 
non-piezoelectric  elastic  half-space  (impedance  Zh ). 


and  piecewise-smooth,  which  implies  that  the  d’Alembert  functions  f\  and  g\  are  also  continuous 
and  piecewise-smooth.  Letting  f[  =  Fx  and  g[  =  G\,  then  the  d’Alembert  solution  9,  together 
with  the  Eq.  7  for  the  stress,  yields 


(Tl(x,t) 


Zl 


F\(t 


x/ci)  -  Gi(t  +  x/ci) 


+  h,D1(x,t )  , 


(10) 


where  z i  =  p i  c ]  is  the  piezo-modified  target  impedance.  Similarly,  the  particle  velocity  v\  in  the 
target  is  given  by 

r\ 

Vi (x,t)  =  =  Fi(t  —  xj Ci)  +  Gi(t  +  x/c1)  .  (11) 

It  is  evident  that  the  stress  and  particle  velocity  in  the  target  are  determined  once  the  functions 
F\ ,  G i,  and  D\  are  determined.  Note  that  discontinuities  in  F\  or  G \  correspond  to  the  passage  of 
shock  waves  at  interior  points  or  the  reflection  of  shock  waves  at  the  boundaries  x  =  0  and  x  —  l. 
This  behavior  is  independent  of  the  back  face  boundary  conditions,  assuming  the  flyer  welds  to 
the  target  at  impact. 
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To  solve  for  F\,  G\,  and  I)  \  requires  consideration  of  the  initial  conditions  (in  this  case 
homogeneous)  and  boundary  conditions  of  the  target.  The  initial  conditions  are 

Mi(x,0)=0,  0)  =  0,  xe[0,Z].  (12) 

Then  Eq.  12i  and  the  stress-strain  Eq.  7  imply  that 

ai(x,  0)  =  Di(x,  0)  =  0  ,  x  G  [0,  /]  .  (13) 


4.  Laplace  Transform  of  the  D’Alembert  Solution 


Solutions  for  the  true  impact  problem  with  appropriate  back  face  boundary  conditions  are 
obtained  in  this  section  by  Laplace  transform  techniques.  The  derivation  of  these  solutions  require 
the  Laplace  transform  of  the  field  equations,  which  is  considered  here.  Let  /  denote  a  real-valued, 
piecewise-continuous  function  defined  at  least  on  [0,  oo).  We  employ  the  right-sided  or  unilateral 
Laplace  transform  C  of  /  defined  as 

OO 

/(»)  =  £{/}(«)  =  £.{/(*)}(»)  =  f  f(t)e-‘dt.  (14) 

0 

Here  s  =  a  +  ice  is  generally  a  complex  number  with  9ft(s)  >  0,  and  the  overbar  on  /  indicates 
that  the  function  is  in  the  Laplace  domain.  The  Laplace  transforms  of  field  variables  such  as 
ui(x,  t ),  <7i(x,  t)  and  vi(x,  t)  are  taken  with  respect  to  the  time  variable  t,  for  example, 

OO 

Vi(x,  s)  —  Ct{(Ji(x,  t)}(s)  —  j  &i(x,  t)  e~stdt  .  (15) 

o 


Applying  Laplace’s  shifting  theorem,  we  see  that 

sx  — SX  — 

Ct{Fi(t  +  x/ci)}  =  eci  -Fi(s)  ,  Ct{Gi(t  -  x/ci)}  =  e  ci  Gi(s)  .  (16) 

Laplace  transformation  of  Eqs.  9,  10,  and  11,  and  using  Eq.  16,  we  obtain  Laplace-domain 
expressions  for  the  displacement,  stress,  and  particle  velocity  in  the  target  in  terms  of  the 
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transforms  of  the  d’Alembert  functions: 


U\  (x,  s)  =  -  e  ci  Fx  (s)  +  e  ci  G1  (s) 

s  L 

r  —sx  —  sx  — 

a i (. x ,  s)  —  Z\  e  ci  F\  (s)  —  e ci  G\ ( s )  |  +  hDi(s)  , 

r  —sx  _  see  _ 

vi(x,  s)  =  e  ci  Fi(s)  +  eci  G*i(s) 


(17) 

(18) 

(19) 


To  find  the  potential  across  the  piezoelectric  target  of  thickness,  /,  integrate  the  electric  field 
E\  =  —  Wj.  Unlike  the  gravitational  field  and  other  field  potentials,  the  minus  sign  is  used  here 
since  the  electric  field  is  “governed  by  Coulomb’s  law  where  elements  of  like  sign  repel,”  cf., 
page  22  of  Sternberg  and  Smith.17 


VM 


l 

j£l(x,t)dx. 

0 


(20) 


Substitution  of  Eq.  2  into  Eq.  20  and  specializing  this  to  a  “short-circuit”  solution,  where  the 
electrical  impedance  between  the  target  impact  face  and  back  face  is  negligible  (see  Fig.  1), 
results  in 

Viit)  =  0  =  —h[ui(l,t)  -  U\ (0,  t)]  +  lDi(t)/e  ,  (21) 


where  the  quantity  in  brackets 
target  surfaces.  Taking  the  Lap 


-  ui(0,t) 


in  Eq.  21  represents  the  relative  motion  of  the 


ace  transform  of  Eq.  21  and  substitution  of  Eq.  17  results  in 


Vi(s)  =  0 


F1(S)(l-e^T)  +  G'1(S)(l 


+  ID i  (s) /e  . 


(22) 


Three  equations  are  required  to  find  the  three  unknown  functions  F\  (s),  Gi(s),  and  I)\  (s).  The 
first  is  given  by  Eq.  21,  the  second  is  given  by  the  boundary  condition  between  the  target  and  the 
half-space  derived  elsewhere,11 

0i (M)  =  zhV!(l,t)  ,  (23) 

and  the  third  equation  is  given  by  a  newly  derived1 1  impact  boundary  condition, 


A)Ul(0,  t)  +  G\  (0,  f )  =  z0V0  . 


(24) 


Although  not  shown  here,  the  elastodynamic  impact  boundary  conditions  (Eqs.  23  and  24)  are 
also  applicable  to  the  piezoelectric  impact  problem  for  a  homogeneous,  elastic,  and  semi-infinite 
flyer  and  half-space  backing  the  target.  Substitution  of  Eqs.  17  through  19  into  the  Laplace 
transforms  of  Eqs.  21,  23,  and  24  permits  unique  determination  of  functions  E\(s),  Gi(s),  and 
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Z)i(s).  Substitution  of  these  functions  into  Eqs.  18,  19,  and  22  allows  determination  of  the  stress, 
velocity,  and  displacement  current  by  using  numerical  ILT  algorithms.  PZT-4  target  mechanical 
and  electrical  properties13  are  given  in  Table  1.  The  elastic  properties  for  the  flyer  and  half-space 
are  identical  and  similar  to  that  of  a  “non-piezoelectric”  steel  material. 


Table  1.  Mechanical  and  electrical  properties. 


Constant 

PZT-413 

Steel 

CE  ( GPa ) 

115.41 

314.89 

zo  =  Poco  {kg  m~2/s 2) 

- 

7850  x  5110 

Zi  =  PI  Cl  {kg  m~2/s2) 

7500  x  4533 

- 

zh  =  z0  {kg  mr2/s2) 

- 

7850  x  5110 

Co  ( farads/m .) 

8.854187817e-12 

- 

er 

663.2 

- 

e  =  e0  er  ( farads/m ) 

5.87209736e-09 

- 

e  ( coul/m 2) 

15.08 

- 

h  =  e/e  {V/m) 

2.5680773e+09 

- 

CD  =  CE  +  eh  {GPa) 

154.136 

- 

Solutions  in  the  time  domain  are  illustrated  in  Figs.  2,  3,  and  4  and  are  found  by  numerical 
inversion  of  the  Laplace  transform  using  the  DAC  algorithm,18  modified  by  Laverty  and 
Gazonas12  to  mitigate  the  effects  of  Gibbs’  phenomenon  using  the  so-called  Lanczos  cr -factors 
with  2048  terms  in  the  Fourier  series  and  a  tolerance  equal  to  10~5.  A  general  compiled 
Mathematica  code  for  the  numerical  inversion  of  functions  with  jump  discontinuities  appears  in 
the  Appendix.  The  stress  in  the  target  in  the  Laplace  domain  s  is  given  by  the  function  containing 
the  exponentials  in  Eq.  A-l.  This  same  code  works  equally  well  for  continuous  functions  and 
users  need  only  replace  the  last  expression  in  Eq.  A-l  with  their  own  Laplace-domain  expression, 
as  well  as  listing  the  constants  appearing  in  their  expression  at  the  beginning  of  the  Module 
command.  Eq.  A- 2  computes  the  numerical  ILT  of  the  function  Stress  and  stores  these  values  in 
the  function  E.  The  Mathematica  compile  function  executes  this  code  in  less  than  8  s  on  a  Dell 
laptop  with  a  64-bit  OS  and  8  GB  of  RAM  running  on  an  Intel(R)  Core(TM)  i5  CPU,  M  560  at 
2.67  GHz.  The  same  code  executes  in  135  s  using  uncompiled  Mathematica  code.  Eq.  A-3  plots 
the  normalized  stress  E  at  the  center  of  the  piezoelectric  target  (see  Fig.  A-l);  the  stress  is 
normalized  by  the  initial  elastic  jump  given  by  the  expression  sjump.  Figs.  2,  3,  and  4  also 
illustrate  solutions  to  the  piezoelectric  impact  problem  using  a  FDTD  method  and  COMSOL 
Multiphysics,  a  commercial  FE  code,13  both  of  which  exhibit  wave  dispersion  effects. 
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Fig.  2.  Normalized  stress  history  at  the  center  of  a  5-mm-thick  piezoelectric 
target  (PZT-4)  backed  by  a  semi-infinite  non-piezoelectric  elastic 
half-space  under  impact  by  a  semi-infinite  non-piezoelectric  elastic  flyer 
with  initial  velocity  Vo  =  5  m/s;  the  normalization  factor  is  the  initial 
elastic  stress  jump  at  impact  given  by  ct(0,  0+)  =  Vo  where  z\e  is 

the  “elastic-only”  impedance  of  the  target,  and 
zie  =  Pi  cie  =  7500  x  3922.8  c.\e  =  \J CE / p\,  see  Table  1. 
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Fig.  3.  Normalized  velocity  history  at  the  center  of  a  5-mm-thick  piezoelectric 
target  (PZT-4)  backed  by  a  semi-infinite  non-piezoelectric  elastic 
half-space  under  impact  by  a  semi-infinite  non-piezoelectric  elastic  flyer 
with  initial  velocity  Vo  =  5  m/s;  the  normalization  factor  is  the  initial 
elastic  velocity  jump  at  impact  given  by  n(0, 0+)  =  ZQ+Z  Vo  ,  where  z\e 
is  the  “elastic-only”  impedance  of  the  target,  and 
zie  =  Pi  cie  =  7500  x  3922.8  c\e  =  \/CE / p\,  see  Table  1. 
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D(t)  coul/m' 


Fig.  4.  Displacement  current  time  history  at  the  center  of  a  5-mm-thick 

piezoelectric  target  (PZT-4)  backed  by  a  semi-infinite  non-piezoelectric 
elastic  half-space  under  impact  by  a  semi-infinite  non-piezoelectric  elastic 
flyer  with  initial  velocity  Vo  =  5  m/s;  the  asymptotic  displacement 
current  Dx{ oo)  =  -  V0  =  -0.0131. 
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5.  Discussion  and  Conclusions 


The  numerical  solutions  reveal  the  nearly  instantaneous  appearance  of  a  tensile  stress  wave  at  the 
center  of  the  target,  which  increases  linearly  in  amplitude  (see  Fig.  2),  and  is  formed  by  the 
£-field  induced  polarization  of  the  material  in  the  initially  unstressed  region  of  the  target.  For 
nonlinear  piezoelectric  media,  Chen  et  al.1  report  that  this  wave  cannot  be  a  shock  wave  but  must 
be  an  acceleration  wave  since  the  wave  speed  depends  upon  a  time-dependent  displacement 
current  D(t)  (see  Chen  et  al.1  equation  (4.12)).  The  nearly  instantaneous  appearance  of  the  stress 
wave  occurs  since  the  speed  of  light  is  about  five  orders  of  magnitude  greater  than  the 
longitudinal  wave  speed  in  the  target.  The  tensile  stress  wave  is  eventually  overcome  by  the 
impact-induced  compressive  stress  wave  at  about  0.5  /as.  We  are  not  aware  of  any  other  “true” 
elastodynamic  impact  solutions  in  the  literature,  which  correctly  treat  the  occurrence  of  multiply 
reflected  transient  waves  in  piezoelectric  targets. 


Numerical  time-domain  solutions  of  stress,  velocity,  and  displacement  current  at  the  center  of  a 
5-mrn -thick  PZT-4  target  are  obtained  via  a  modified-DAC  algorithm,  and  compare  well  with 
those  obtained  using  both  a  FDTD  method  and  the  commercial  FE  code  COMSOL 
Multiphysics.13  Although  not  shown  in  this  report,  the  final  value  theorem  (FVT)  can  be  used  to 
verify  that  the  asymptotic  (long-time)  stress  and  velocity  in  the  piezoelectric  target  are 
independent  of  the  target’s  elastic  or  piezoelectric  properties  and  only  the  elastic  properties  of  the 
flyer  and  half-space,  a  result  also  proven  using  the  time-domain  solutions  obtained  for  purely 
elastic  targets.11  The  FVT  is  also  used  to  determine  a  closed-form  expression  for  the  asymptotic 
displacement  current, 


-Di(oo)  = 

illustrated  in  Fig.  4,  where  ay  (a;, 


_ Zp  zh  he  _  o~i(g,oo) 

(z0  +  zh)  (ci  zi  -  h2e)  0  (h  -  ^-) 

oo)  is  the  asymptotic  target  stress.11 


(25) 
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Appendix.  Mathematica  Code  for  Determining  the  Numerical  ILT  Using  the 
Modified-DAC  Algorithm 


The  following  is  the  Mathematica19  code  for  determining  the  numerical  inverse  Laplace  transform 
(ILT)  using  the  modified  Dubner-Abate-Crump  (DAC)  algorithm.12  Figure  A-l  shows  the 
normalized  stress  time  history  of  the  target. 


Stress  =  Compile[{{s,  jComplex}},  ModuIe[{/,  x,  eo,  Ce,  ers,  e2,  e,  h,  M,  p,  cl,  VO, 
zO,  zl,  zh,  a,  b,c},l  =  5  *  10“3;  x  =  1/ 2;  eo  =  8.854187817  *  10~12;  Ce  =  115.41  *  109; 
ers  =  663.2;  e2  =  15.08;  e  =  eo  *  ers ;  h  =  e2/e;  M  —  Ce  +  e 2  *  h;  p  =  7500; 


cl  =  J  — ;  VO  =  5;  zO  =  7850  *  5110;  zl  =  p  *  cl;  zh  =  zO; 
(z0  —  zl)  (zl  —  zh)  zO  —  2zl  +  zh 

a=W^T)j^rb=z0+2zl+zh’c= — a — ; 

/  a(t+a:)  a(l  +  2x)  is  »\ 

V0z0e~zi  f  [zz*-e  cl  +  zle  a  +  zlea  —  [zl  —  zh)ea  J 
s  h  (-e^  +  Islsk  +  c)  +  ^+zi)(zi+zh)ls  ^  + 

'-e^(zl  +  zh)  -  zle^  (l  +  Is&p f1)  +  zle% 

b  +  4 zlM  +  +  (zO+zl)(z^l+zh)  ls  ^  + 


I  42leST 


]]  (A-l) 
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E  =  Compile[{},  Module [ { sumcp ,  sumsp,  tf,  tol ,  a,  tot,  cO,  r,  dt,  a,  arg,  coeffs, 
cospart,  cl,  c2,  sinpart,  kpd,  t,  kpt};  tf  =  5  *  1CT6;  tol  =  0.00001; 

a  =  —  Log  [toZ]  _  ^  _  2Q4g.  cq  _  stress  [a];  r  =  1000;  dt  =  — ; 

2  tj  t 

kn  kn 

Table [sumcp  =  0.;  sumsp  =  0.;  Do  [kpd  =  — ;  kpt  =  — ; 

o’  =  '  -  7  ^ ;  coef  fs  =  Stressfa  +  i  *  kpt];  cl  =  Re[coe//sl;  c2  =  Im[coe//sl; 

kpd 

cospart  =  (TclCosft  *  kpt] ;  sumcp  =  cospart  +  sumcp; 

sinpart  =  <7c2Sin[t  *  /cpt]);  sumsp  =  sinpart  +  sumsp,  {k,  1,  tot}]; 

Exp  [a  *  t]  ( v  +  sumcp  —  sumsp )  r 
— - -  2  tf - - - — ,  {t,  dt,  tf,  dt}]], 

CompilationOptions  — *  {" Ini ineExtcrnal Definitions"  — *  True, 
"InlineCompiledFunctions“  — >•  True}] ; 


VO  =  5;  p  =  7500;  zO  =  p  *  5110;  Ce  =  115.41  *  109;  zn  =  sj p  *  Ce; 

zOzn 


sjump  =  V0 - 


-;  ListPlot[E[  ]/ sjump,  Joined  — >•  True, 


zO  +  zn 

BaseStyle  — >  (FontFamily  — y  "Times",  FontSize  — >  12}, 
PlotRange  — >  Full,  DataRange  — y  (0,  5  *  10~6}, 

.  cr(Z/2 ,t) 


AxesFabel  — >  ("Time  (s)", 


o-(0,0) 


‘}] 


(A-2) 


(A-3) 
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o' (1/2,  t) 


cr  (0,  0) 


Tjme  (s) 


Fig.  A-l.  Normalized  stress  time  history  at  the  center  of  a  5-mm-thick 

piezoelectric  target  (PZT-4)  backed  by  a  semi-infinite  non-piezoelectric 
elastic  half-space  under  impact  by  a  semi-infinite  non-piezoelectric 
elastic  flyer  with  initial  velocity  Vq  =  5  m/s. 
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